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The large N commensurate dirty boson model, in both the weakly and strongly 
commensurate cases, is considered via a perturbative renormalization group treatment. 
In the weakly commensurate case, there exists a fixed line under RG flow, with varying 
amounts of disorder along the line. Including 1/N corrections causes the system to flow 
to strong disorder, indicating that the model does not have a phase transition pertur- 
batively connected to the Mott Insulator-Superfluid (MI-SF) transition. I discuss the 
qualitative effects of instantons on the low energy density of excitations. In the strongly 
commensurate case, a fixed point found previously is considered and results are obtained 
for higher moments of the correlation functions. To lowest order, correlation functions 
have a log-normal distribution. Finally, I prove two interesting theorems for large N 
vector models with disorder, relevant to the problem of replica symmetry breaking and 
frustration in such systems. 

I. INTRODUCTION 

The dirty boson problem, a problem of repulsively interacting bosons in a random potential, has 
been the subject of much theoretical work In the zero temperature quantum problem, the system 
can undergo a phase transition between a Bose glass phase and a superfluid phase. An action that may 
be used to describe this transition is 

d d x dt (+d x lj>(x, t)d x (j>{x, t) + d t 0{x, t)d t <j>(x, t) + w(x)&(x, t)d t <t>(x, t) ~ U{x]4>(f> + g(4>cj>) 2 ^ (1) 

This action describes a system with a number of phases. In the pure case, in which U(x) is 
a constant, and w(x) vanishes everywhere (the commensurate case), there is a phase transition from a 
gapped Mott insulator phase to a superfluid phase as U increases. If we consider a case in which U(x) is 
not constant, but w(x) still vanishes everywhere, there will be a sequence of transitions from a gapped 
Mott insulator, to a gapless (but with exponentially vanishing low energy density of states) Griffiths 
phase, and then to a superfluid phase. The Griffiths phase occurs due to the possibility of large, rare 
regions in which fluctuations in U(x) make the system locally appear superfluid. 

If we consider the case in which both U(x) and w(x) are fluctuating, there will also appear a 
Bose glass phase in which there is a density of states tending to a constant at low energy and an infinite 
superfluid susceptibility ||. The physical basis for this phase is the existence of localized states, in which 
a competition between chemical potential and repulsion causes the system to desire a certain number of 
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particles to occupy each localized state. There exist excitations involving adding or removing one particle 
from these states, and these excitations lead to the diverging susceptibility. However, it is clear that the 
Bose glass phase is very similar to the Griffiths phase, in that both involve regions of finite size, fn the 
Griffiths phase one needs regions of arbitrarily large size, while in the Bose glass phase one only needs 
regions large enough to support a localized state, with a nonvanishing number of particles occupying that 
state, in order to produce the diverging susceptibility characteristic of the phase. In a system in which 
the disorder is irrelevant at the pure fixed point, so that the fluctuations in U(x) and w(x) scale to zero, 
one will still find a Bose glass phase as there will, with low probability, exist regions that can give rise 
to these localized states. Thus, the interesting question to answer is not whether the Bose glass phase 
exists, but whether there exists a fixed point at which fluctuations in w(x) are weak, so that the critical 
exponents are near those of the pure Ml-SF transition. The most likely alternative would be governed 
by the scaling theory of Fisher et. al., which has very different critical exponents We will refer to 

this scaling theory as the phase-only transition, as one assumes fluctuations in the amplitude of the order 
parameter are irrelevant at the critical point. 

Recently, a large N generalization of equation (|]) was considered in the restricted case w(x) = 
j|. We will refer to this case as the strongly commensurate case, while the situation in which w(x) 
vanishes on average, but has nonvanishing fluctuations will be known as the weakly commensurate case. 
We consider a system defined by the partition function 

/ (n s @i( x > wo*. *) - N °*{*)j) i#~ s ( 2 ) 

^ x,t 

where 

S= [ d d xdt(d x J l (x,t)d x ^\x,t) + dtMx,t)d t ^\x,t) + w(x)^ l (x,t)d t <p l {x,t)) (3) 



Here, we have, for technical simplicity later, replaced the quartic interaction by a ^-function. For most 
of the paper, a J-function interaction will be used. However, for generality in the last section, we will 
return to quartic interactions. The disorder in U(x) will be replaced, in the 5-function case, by weak 
fluctuations in a 2 . We will consider a 2 = a 2 + Sa 2 where a\ is a constant piece used to drive the system 
through the phase transition and Sa 2 is a fluctuating piece. 

The advantage of the large N formulation of the problem is that for any fixed realization of the 
problem one may exactly solve by system by finding the solution of the self-consistency equation 

a 2 (x) = (x,t = 0\(-d x r -d 2 + w(x)d t + \{x))- 1 \x,t = 0) (4) 



or 



<t 2 (x) = / du{x,t = 0\{-dl+oj 2 +iw(x)uj + \{x))- 1 \x,t = 0) (5) 



where A(x) is a Lagrange multiplier field for enforce the 5-function constraint on the length of the spins. 
After solving the self-consistency equation, any correlation function can be found simply by finding the 
Green's function of a non-interacting field with action J d d x dt 4>(—d 2 — d 2 + w(x)d t + A(x))(f>. 
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In equation (Q), we assume that the Green's function on the right-hand side has been renormalized 
by subtracting a divergent quantity. Specifically, we will take a Pauli-Villars regularization for the Green's 
function, and take the regulator mass to be very large, while adding an appropriate divergent constant to 
a 2 on the left-hand side. The cutoff for the regulator is completely different from the cutoff for fluctuations 
in da 2 that will be used for the RG later; the cutoff for the regulator will be much larger than the cutoff 
for fluctuations in Sa 2 and w(x), and will be unchanged under the RG. 

In the previous work, exact results were obtained for the critical exponents for average quantities. 
In the present paper, I will first consider the problem in which w{x) ^ 0, although w(x) vanishes on 
average (the weakly commensurate dirty boson problem). I will consider general problems of the large 
N system with terms linear in the time derivative. Next, a lowest order perturbative RG treatment will 
be used to consider critical behavior. Instanton corrections to the perturbative treatment will be briefly 
discussed after that. Returning to the strongly commensurate case, previous results on the fixed point 
will be extended to give results for higher moments of the correlation functions. Finally, as a technical 
aside, we consider the large N self-consistency equation in frustrated systems, and demonstrate that 
the self-consistency equation always has a unique solution, as well as considering the number of spin 
components needed to form a classical ground state in frustrated systems. 



II. BOSE GLASS IN THE LARGE TV LIMIT 



Consider the following simple 0+1 dimensional problem, at zero temperature (/3 — > oo): 

^(*)]^*(^(*)^(*)-^^ { "*^ +> ***' },,, (6) 



The solution of this problem in the large N limit via a self-consistency equation requires finding a A such 
that 

^ I (7) 

27T iAiO + LU 2 + A 

Then, by contour integration, for A > 0, we find a 2 = V J +4X ■ Then, 

Unfortunately, this result for A leads to A becoming negative for sufficiently large A. Although pertur- 
bation theory will not see a problem, this is the signal for the Bose glass phase. One must separate the 
self-consistency equation into two parts, one containing an integral over non-zero uj, and one containing 
the term for zero u>. One finds (for finite (3) 

2 _ f duj 1 11 

a ~ J ~2tt iAu + + A + A ( ' 
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Then, the self-consistency equation can always be solved using positive A, but in the zero-temperature 
limit of the problem one will find that one needs A to be of order -g, and at zero- temperature, there will 
appear a zero energy state. 

Considering the original statistical mechanics problem of equations (|l]j|), one will expect to see 
some non-zero density of these zero energy states, indicating the presence of a gapless Bose glass phase, 
with diverging superfluid susceptibility. Even if the fluctuations in w(x) vanish at the fixed point, when 
the critical exponents are unchanged from the strongly commensurate problem, such zero energy states 
will exist as Griffiths effects leading to the appearance of a Bose glass phase near the superfluid transition. 

To perform a renormalization group treatment of the model, we will first proceed in a perturbative 
fashion in the next section, ignoring such zero energy effects. For small fluctuations in w(x), these zero 
energy states will be exponentially suppressed as will be discussed in the section after that. 



III. PERTURBATIVE RG 

We will follow the RG techniques used in previous work on the large N problem PJ. We will 
work near 2 + 1 dimensions, specifically we will have d = 2 + e space dimensions and 1 time dimension. 
We will at first work to one-loop in perturbation theory, which will, in some cases, correspond to lowest 
order in e. Some results will be extended to all orders. A fixed line is found for the large N system. The 
fixed points are destroyed by 1/N corrections. 

The RG is defined as follows: start with a system containing fluctuations in Sa 2 and w(x) up 
to some wavevector A. Remove the high wavevector fluctuations in these 6a 2 and w(x) to obtain a new 
system, with renormalized gradient terms d 2 and d 2 in the action, as well as renormalized low wavevector 
Sa 2 and w{x) terms. Do this procedure so as to preserve the average low momentum Green's function, 
as well as the low wavevector fluctuations in \{x). See previous work [Q for more details. 

If we are working in 2 + e space dimensions, and 1 time dimension, we can easily work out 
the naive scaling dimensions of the disorder. One finds that if we assume Gaussian fluctuations in the 
disorder, with 

(8a 2 {p)5a 2 {q)) = (2n) 2 S(p + q)S (10) 

(w(p)w(q)) = (2tt) 2 S( P + q)W (11) 

then S scales as length to the power d — 2 = e, while W scales as length to the power 2 — d = — e. So, for 
d > 2 we find that disorder in a 2 is relevant at the pure fixed point, while disorder in w(x) is irrelevant. 
For d < 2 this is reversed. 

Previously, a lowest order in e calculation Q considering only disorder in 6a 2 gave the following 
results. For a given problem, with fluctuations in Sa 2 up to wavevector A, and self-consistency equation 
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a 2 (x) = j du{x,t = Q\(-dl + u 2 + X{x)^j 1 )\x,t = Q) (12) 



one could define another problem, with fluctuations in Sa 2 only up to A — SA, with self-consistency 
equation 

(1 - S -±c 3 L)a 2 (x) = J S D - d ^ {x,t= 0| ( (1 + 6 -±c 2 L)d 2 x + (1 + 6 -^c 3 L)(uj 2 + A(aO)) _1 )M = 0} 

(13) 

such that the Green's function computed from the second self-consistency equation agrees with the Green's 
function computed from the first self-consistency equation averaged over disorder at large wavevector. 
Here we define £ = c 2 A s ~ 2D S where c\ = ^°^ 2 r ( 2 -D/2)r^^D/2-i) ' C2 = (l~4/d)c 3 , andc 3 = 2 ^^ A d ~ 4 . 

The results above, to one loop, were obtained by considering the large wavevector fluctuations in 
A due to the large wavevector fluctuations in Sa 2 , and then finding how they renormalize the self-energy 
and vertex. To lowest order, one obtains the fluctuations in A by inverting a polarization bubble. That 
is, one expands the self-consistency equation to linear order in A to solve for large wavevector fluctuations 
in A as a function of fluctuations in Sa 2 . One finds then that 

6<j 2 (p) = c~ 1 p D - i \(p) + ... (14) 

From this, we obtain fluctuations in A at wavevector A which, to lowest order, are Gaussian with mean- 
square L. See figures 1,2, and 3. For more details, see previous work p|. 

It may easily be seen that, to lowest order, the addition of the term w(x) does not produce any 
additional large wavevector fluctuations in A, as equation (|lj) is still true to lowest order in w(x) and 
X(x). However, the term w(x) can produce a renormalization of the self-energy. See figure 4. The result 
is to produce a term in the self-energy equal to 



= -5Auj 2 J ^ k —j— W (15) 

This is equal to 

6A -u 2 Wc A + ... (16) 



' (p + kf 

fc 2 =A 2 



A 

where c 4 = A d ^ 2 2 ^J 2y 

There is one other term that must be included in the RG flow equations at this order. The 
fluctuations in A due to the fluctuations in Sa 2 can renormalize the vertex involving w(x). See figure 
5. This will change the term w{x)dt in the self-consistency equation to w(x)dt{l + ^-c 3 L). Note that 
the renormalization of the w(x) term is equal to the renormalization of the to 2 and X(x) terms in the 
self-consistency equation. 

Putting all the terms together, we find that with a lowered cutoff A — SA, the renormalized theory 
is described by the new self-consistency equation 
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(l-5 3 )a 2 (x) = J du (x,t = 0| ( (1 + S 2 )d 2 x + (1 + 6 4 )lj 2 + i(l + S 3 )w(x)u + (1 + 5 3 )\(x)) 1 )\x,t = 0) 

(17) 

where S 3 = ^-c 3 L, 5 2 = ^-c 2 L, and £ 4 = <5 3 + ^-Wc^. Rescaling oj by (1 + &2 ~ Si ) to make the coefficients 
in front of the oj 2 and d 2 terms the same, rescaling A, and then rescaling the spatial scale to return the 
cutoff to A we find 

»•(.) =/*,<,. t- 0| (-* + „» + + AW) - ')|x. . = 0) ,18) 



where 



a 2 = (l-5 3 + 5 2 + -^ + (l + e) T ) f 7 2 (x) (19) 



*(x) =»(! + 6 3 - 6 2 + 5jL ^ ± )w(x) (20) 



From this, we extract RG flow equations for Oq, S 1 , and W. The result is 

^ mCr r t tj/ C4 



rflnA 

dlnT^ 
rflnA 



l + e-c 3 L + c 2 L + Wj (21) 
= e - 2c 3 L + 2c 2 L + Wc 4 (22) 
= -e + 2c 3 L - 2c 2 L - Wc 4 (23) 



The renormalization group flow has a fixed line, as the product SW is invariant under the RG 
flow. It may be verified that the ratio S/W has a stable fixed point under RG flow for any e and any 
value of SW. Further, it may be seen that the critical exponent v on the fixed line is given by vd = 2, as 
if S is constant under RG flow, then ctq has = 1 + f d/2. Later, we will consider the effect of 1/N 

corrections. 

First, note that the line is peculiar to having one time dimension. For fewer than one time 
dimension, there will be a stable fixed point at W — 0, which is attractive in the W direction. Thus, 
in the framework of a double-dimensional expansion, one may not sec problems at low orders, as the 
fixed point has nice behavior for small numbers of time dimensions. Compare to results in the double- 
dimensional expansion ||. 

Further, the presence of the fixed line only required that the renormalization of the w{x)dt vertex 
was equal to the renormalization of the vertex on the left-hand side of the self-consistency equation 
defining a 2 . This equality will persist to all orders in a loopwise expansion via a Ward identity. Thus, 
we expect that the fixed line is an exact property of the large N theory. 
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Let us consider the effect of 1/N corrections on this line. To lowest order in 1/N, for weak 
disorder, the 1/7V corrections only modify the naive scaling dimensions in the RG flow. The scaling 
dimension of 4>dt4> is not changed under 1/N corrections. However, the scaling dimension of c/xf) is 
changed by an amount r\ = -Jp-^f- Thus, 1/N corrections will change the RG equations to 

dlnS" 



dlnW 



= 2ri + e- 2c 3 L + 2c 2 L + Wca (24) 
= -e + 2c 3 L - 2c 2 L - Wc 4 (25) 



rflnA 

Then, we find that SW is growing under the RG flow, and the system goes off to a different fixed point. 
The most reasonable guess then is that in a system with finite N (including physical systems with N = 1 ), 
the transition is not near the MI-SF transition, but is instead of another type, perhaps the phase-only 
transition. Other authors have shown that, in some cases, the phase-only transition is stable against 
weak commensuration effects 101 . 



IV. INSTANTON CALCULATIONS 



Unfortunately, the ability to carry out instanton calculations in this system is rather limited. 
It will not be possible to calculate the action for the instanton with any precision, but we will at least 
present some arguments about the behavior of the instanton. The idea of the calculation is to look for 
configuration of w(x) and <J 2 (x) (these configurations are the "instantons"), such that the self-consistency 
equation cannot be solved without including contributions from zero energy states, as discussed in section 
2. Let us first consider the case in 2 + 1 dimensions. Let us assume that we try to produce such states 
in a region of linear size L. Looking at the lowest energy state in this region, one would expect that 
the contribution of spatial gradient terms in the action would lead to an energy scale of order L^ 1 . The 
linear term w(x)u in the action will become important, and produce such a zero energy state, when w(x)lo 
becomes of order uj 2 . This implies occurs when w(x) is of order u, which implies w(x) ~ L^ 1 . For some 
appropriate configuration of w(x), assuming quadratic fluctuations in w(x) with strength of order W, 
we will have an action Smstantcm oc J - ' ffl d 2 x. Thus, these configurations will occur with exponentially 
small probability e ~ ,Sinatanton for weak disorder in w(x). 

Away from 2 + 1 dimensions, one will find that the action for the instanton, ignoring fluctuation 
corrections, is dependent on scale. For d > 2 it is increasing as the scale increases, indicating that 
large instantons are not present. For d < 2, it is decreasing as the scale increases, indicating that large 
instantons are easy to produce. This is simply a way of restating the fact that fluctuations in w(x) are, at 
the pure fixed point, irrelevant for d > 2 and relevant for d < 2. It is to be expected that corrections due 
to fluctuations as considered in the renormalization group of the previous section will make the action 
for the instanton scale invariant. However, since we do not fully understand how to calculate instanton 
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corrections even in the simplest d — 2 case, the task of combining instanton and fluctuation corrections 
is presently hopeless. 



V. HIGHER MOMENTS OF THE GREEN'S FUNCTION 

Having considered the weakly commensurate case, and found no fixed point in physical systems, 
we return to the strongly commensurate case with W = 0, and consider the behavior of higher moments 
of the Green's function. A lowest order calculation will show log-normal fluctuations in the Green's 
function. 

Let us first consider the second moment of the Green's function. That is, we would like to compute 
the disorder average of the square of the Green's function between two points, which we may write as 
(G(0,x) 2 }. We may Fourier transform the square to obtain (G 2 (p,oj)}. Now, one may, when averaging 
over disorder, include terms in which disorder averages connect the two separate Green's function. At 
lowest order, there is no low-momentum renormalization of the two Green's function propagator, beyond 
that due to the renormalization of each Green's function separately. See figure 6. That is, if one imagines 
the two Green's functions entering some diagram, with both Green's functions at low momentum, going 
through a sequence of scatterings, and exiting, again with both Green's functions at low momentum, 
one does not, to lowest order, find any contribution with lines connecting the two Green's function. The 
reason for this is that at this order we will only join the Green's functions with a single line, along which 
one must have momentum transfer of order A. This then requires that some of the ingoing or outgoing 
momenta must be of order A. 

One does, however, find a contribution which we may call a renormalization of the vertex. See 
figure 7. In order to find the second moment of the Green's function, one must start both Green's 
functions at one point, and end both Green's functions at another point. Near the point at which both 
Green's functions start, one may connect both lines with a single scattering off of A, at high wavevector of 
order A. Then, one can have a large momentum of order A circulating around the loop formed, while the 
two lines that leave to connect to the rest of the diagram still have low momentum. This then replaces 
the two Green's function vertex, which we will refer to as V2, by a renormalized vertex. 

The result of the above contribution is that the two Green's function vertex V2 is renormalized 
under RG flow as 

^ = C3i = c,c ? A— 8 (26) 

Then, the second moment of the Green's function, at momentum scale p is given in terms of the first 
moment by 

(G 2 (p,c))cx(G(p,c))V 2c3i (27) 
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Note the factor of 2 in front of C3L, as the second moment of the Green's function gets renormalized at 
both vertices. One must insert the value of L at the fixed point into the above equation to obtain the 
behavior of the second moment. 

For higher moments, the calculation is similar. In this case, one must, for the n-th moment, 
renormalize a vertex V n . The result is 

rfln ^ "(» - 1) T n{n-l) 2 S _ 2D 

dmA = * L = ~^r~ C3ClA S (28) 

The factor arises as at each stage of the RG one may connect any one of the lines in the vertex 

V n to any other line in the vertex. There are "^" 2 ~ 1 ^ ways to do this. Then one finds 

(G"(p, u)) oc (G{p, u)) n p- n(n - 1)c * L (29) 

This result for the behavior of the higher moments of the Green's function is quite typical 
for disordered systems. Compare for example to the results on 2-dimensional Potts models |(|. From 
the results for the moments of the Green's function one may, under mild assumptions, determine the 
distribution function of the Green's function. This distribution function is the probability that, for a 
given realization of disorder, the Green's function between two points assumes a specific value. From 
the result for the moments given above one finds that the distribution function is log-normal. That is, 
the log of the function has Gaussian fluctuations. Physically this should be expected from any lowest 
order calculation, as lowest order calculations generally treat momentum scales hierarchically, and one 
is simply finding that at each scale there are random multiplicative corrections to the Green's function, 
causing the log of the Green's function to obey a random walk as length scale is increased. 



VI. GLASSY BEHAVIOR IN THE LARGE TV LIMIT 



First, we would like to demonstrate that, in the large N limit, the self-consistency equation 
always has a unique solution. For generality, we consider here the case of quartic interactions instead 
of (5-function interactions. In the absence of terms linear in u>, uniqueness is clear on physical grounds, 
for the models considered above in which the coupling between neighboring fields (f> is ferromagnetic and 
unfrustrated. However, we will show this to be true for any coupling between neighboring fields and in 
the presence of terms linear in to. 

Note that, for finite N, the terms linear in oj lead to frustration. Consider an N = 1 system 
with a finite number of sites. Assume that there is no hopping between sites, but there is some repulsion 
between sites due to a quartic term. Let there be terms linear in u> in the action, but no terms quadratic 
in oj. The states of the theory are then determined by how many particles occupy each site. The repulsion 
leads to an effective anti-ferromagnetic interaction, in the case in which each site has zero or one particles 
and we imagine one particle to represent spin up and no particles to represent spin down. This can then 
produce frustration. Compare to the Coulomb gap problem in localized electron systems pj . 
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However, physically speaking, as N — > oo, the discreteness of particle number on each site 
disappears, and the system becomes unfrustrated. We will now show this precisely. 

Consider a problem at non-zero temperature, so that there is a sum over frequencies u>. Consider 
an arbitrary single particle Hamiltonian Hq, defined on a V-site lattice, so that the self-consistency 
equation involves finding Xi, where i ranges from 1 to V, such that 

o 2 + M U*i = J2(i\( H o + ^ + iA.lu + BiW 2 )- 1 !*) (30) 

3 u 

Here, o 2 is a function of site i, and A\ and are functions of site i defining the local value of the linear 
and quadratic terms in the frequency. The matrix Mij is included to represent the effects of a quartic 
interaction. For the problem to be physically well defined, M+j must be positive definite. 

The proof that equation (|3(]) has only one solution proceeds in two steps. First we note that if 
Hq vanishes, then the equation obviously only has one solution with A > 0. Next, we will show that as 
Hq varies, Aj varies smoothly, and therefore any arbitrary Hq can be deformed smoothly into a vanishing 
Ho, leading to a unique solution for Xi for arbitrary Hq- 

Consider small changes 5Hq and SXi. In order for the self-consistency equation to remain true, 
if we define 

iH = -^2(i\(H + Xi + iAiU + BiU) 2 )- l 8H Q {Ho + A, + iA lL u + BiU) 2 )- l \i) (31) 
we must have 

«i = ^2(i\( H o + A,: + iAiLo + BiUpy-SXiiHo + A, + iAiUJ + BiuP)- l \i) + ^ MijXj (32) 

w 3 

The right hand side of equation ( |32] ) defines a linear function on SXi. If it can be shown that this function 
is invertible, then the theorem will follow. However, we have that 

Tr((5A; J2{H + X t + iA t uj + Biu; 2 )- 1 SXi{H + A. t + iA.uj + BiU 2 )- 1 ) > (33) 

OJ 

due to the well known fact that second order perturbation theory always reduces the free energy of 
a quantum mechanical system with a Hermitian Hamiltonian at finite temperature. We also have, as 
discussed above, that is positive definite. Thus, the linear function on SXi defined above is a sum of 
positive definite functions, and hence positive definite. Therefore, it is invertible and the desired result 
follows. 

This result is interesting considering the phenomenon of replica symmetry breaking. It has been 
noticed by several authors that large N infinite-range spin glass models do not exhibit replica symmetry 
breaking within a meanfield approximation, both in the classical and quantum cases § . Although those 
calculations were based on the absence of unstable directions, in the large N limit, for fluctuations about 
the replica symmetric state, it is possible that the real reason for the absence of replica symmetry breaking 
is the uniqueness of the solution of the self-consistency equation, as shown above. 
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A second interesting question, having begun to consider possible glassy behavior in the large N 
limit, has to do with the nature of the ground state in the classical limit. If we drop all terms in u), to 
produce a classical problem, and ask for the classical ground state, for some arbitrary bare Hamiltonian 
Ho, one may ask how many of the N available spin components will be used. 

In this case, consider Hamiltonian Ho, which is a V^-by-V^ matrix in the case where there are V 
sites. First consider the case in which H is a real Hcrmitian matrix. Since we are considering arbitrary 
Hamiltonians Ho, we can, without loss of generality, constrain all spins to be the same length. We can 
find the classical ground state by looking for solutions of the self-consistency equation 

<7 2 = (i|(flb + A i )- 1 |i) (34) 

in the limit as a 2 — > oo. 

In this limit, the right-hand side will be dominated by zero energy states (more precisely, states 
that tend to zero energy as a 2 tends to infinity) of the operator Hq + A. If the system has k of these 
states, the ground state of the system will use k of the spin components. If the system needs to use all 
k of these components to form a ground state, that is, ignoring the case in which a state using k spin 
components is degenerate with a state using fewer components, then even under small deformations of 
Hq the system will use k spin components in the ground state. Then, under these small deformations, 
Hq + A will still have k zero eigenvalues. To produce k zero eigenvalues for all real Hermitian matrices in 
a neighborhood of a given Hermitian matrix Ho requires k(k + l)/2 free parameters. The elements of A 
provide these parameters. Since there V of these elements, we find that k(k + l)/2 < V, and the number 
of spin components needed to form the classical ground state is at most \/2V. 

If Hq were an arbitrary Hermitian matrix, with complex elements, or a symplectic matrix, one 
would find a similar result, with k still at most order y/V, although the factor 2 would change. This is 
analogous to the different universality classes in random matrix theory Q . Finally, we make one note on 
the number of parameters available to solve the self-consistency equation. There are V free parameters. 
However, self-consistency requires solving V independent equations, so the number of variables matches 
the number of equations. By considering the number of parameters requires to produce zero eigenvalues 
of H + A, we were able to obtain a bound on the number of zero eigenvalues. Still, one might wonder if 
there are enough free parameters to produce multiple zero eigenvalues and still solve the self-consistency 
equations, as it appears that one would then need k{k + l)/2 + V free parameters. However, if there 
are k zero eigenvalues, by considering the different ways of populating the zero energy states (that is, 
considering the different ways in which the eigenvalues tend towards zero as a 2 tends toward infinity) 
one obtains an additional k{k + l)/2 parameters, so the number of parameters available always matches 
the number of equations. 

We can extend this theorem to look at metastable states. Suppose a configuration of spins is 
a local extremum of the energy Ho, for fixed length of spins. Then, since the derivative of the energy 
vanishes, one finds that a matrix (Ho + Xi) must have a number of zero eigenvalues equal to the number of 
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spin components used. Suppose that for small deformations of H there is still a nearby local minimum, 
as one would like to require for a stable state. Then, one can argue that the number of spin components 
k used in the state obeys k(k + l)/2 < V. 

This second theorem may be of interested in considering the onset of replica symmetry breaking. 
If we have a system in a large volume and large N limit, one must ask in which order the limits are taken. 
If the N — > oo limit is taken first, there will be no replica symmetry breaking. However, if the infinite 
volume limit is taken first, there may be replica symmetry breaking. If one has N > 2fc max , where fc max 
is the largest k such that k(k + l)/2 < V, then there arc no local minima other than the ground state. 
This follows because, as shown above, a local extremum of the energy, fa, will use at most k max spin 
components, and the ground state, 4>f ', can be constructed using a different set of fc max spin compoents. 
Then starting from fa, one finds that deforming the state along the path \J\ — 8 2 fa + 5<j)f as S goes from 
to 1 provides an unstable direction for fluctuations. 

VII. CONCLUSION 

In conclusion, we have considered the large N dirty boson model, including the effects of local 
incommensuration (the terms linear in lu). In the large N limit, a fixed line under RG is found, but 
is destabilized by including 1/iV corrections. This suggests that the phase transition in experimental 
(N — 1) systems is of the phase-only type, instead of the MI-SF type. 

There is a problem with local incommensuration in a perturbative approach, as discussed in the 
section on the Bose glass in the large TV limit and the section on instanton calculations. One would like 
a quantitative method of assessing the results of the instantons, although this is largely a technical issue, 
as it appears that there are no accessible fixed points in the RG using this approach. 

In the strongly commensurate case, it has been shown that one can calculate higher moments of 
the correlation functions. The result shows that the correlation functions have a log-normal distribution. 

The two theorems proved in the last section give useful information on the relevance of the large 
N expansion in frustrated problems. It would be interesting to use these results as a starting point for a 
better understanding of replica symmetry breaking. 

The large N approximation has been a useful approximation for pure, unfrustrated systems. It 
is hoped that it may become as useful for disordered interacting systems. 
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FIG. 1. Polarization bubble. Thick lines represent either scattering vertex off A or scattering vertex used to 
define a 2 in self-consistency equation. 




FIG. 2. Self-energy correction due to fluctuations in A. Joining the thick lines in a loop denotes averaging A 
over disorder in a 2 . Momentum of order A flows around loop. 




FIG. 3. Vertex correction due to fluctuations in A. This represents both renormalization of vertex defining 
scattering off of A and renormalization of vertex defining a 2 in self-consistency equation. 
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FIG. 4. Renormalization of uj 2 term due to fluctuations in w(x). The vertices with thin lines represent scattering 
off of w(x)dt, while the joining of the vertices represents averaging w(x) over disorder. 



FIG. 5. Renormalization of vertex for w(x) due to fluctuations in A. 



FIG. 6. Possible contribution to propagation of two Green's functions, both in same realization of disorder. This 
is used to compute second moment of the Green's function. This diagram does not lead to any renormalization 
of low momentum behavior. 
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FIG. 7. Renormalization of vertex in computing higher moments of Green's function. Two Green's functions 
start at the same point. After Fourier transforming, this implies that they start with given total momentum. By 
including fluctuations in A, with momentum of order A running around the loop, one can define a renormalized 
vertex. 
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